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Abstract 

The present work studies the Ghatak-Sherrington (GS) model with the presence of a longi- 
tudinal magnetic random field (RF) /i, following a bimodal distribution. The model considers 
a random bond interaction Ji j which follows a Gaussian distribution with mean JqIN and vari- 
ance IN. This allows us to introduce the bond disorder strength parameter 7/70 to probe the 
combined effects of disorder coming from the random bond and the discrete RF over unusual 
phase transitions known as inverse transitions (ITs). The results within a mean field approxi- 
mation indicate that these two types of disorder have complete distinct roles for the ITs. They 
indicate that bond disorder creates the necessary conditions for the presence of inverse freezing 
or even inverse melting depending on the bond disorder strength, while the RF tends to enforce 
mechanisms that destroy the ITs. 
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1. Introduction 

Inverse Transitions (ITs) are counterintuitive reversible transitions, where the phase usually 
considered the more ordered has more entropy than the disordered one. ITs have been observed 
experimentally in very differentphysical systems, ranging from gold nanoclusters to high-Tc 
superconductors ||jl|2l|3l|4l|5l|aLZl|8l|9|]. Such diversity of physical systems displaying ITs have 
also triggered intense theoretical interest in which several models have been studied, in particular, 
magnetic ones because of its simplicity 

lis [U E [il [3 [B E I ITLlISr . From those magnetic 



models, it is suggested that disorder and frustration would be key ingredients to produce ITs, 
freezing (IF) or even melting (IM) {l^. However, most of the previously mentioned studies 
using magnetic models suffer from a limitation. The disorder comes basically from random 
bonds (RB) between spins. Therefore, it is an open issue whether or not other forms of disorder 
can contribute to produce ITs. 
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The family of Blume-Emery-Griffiths model (l^ as the Blume-Capel (EC) ^ and Ghatak- 
Sherrington (GS) ones 112 ill have been quite useful to study the ITs. Therefore, they can be 
also helpful to answer the issue mentioned previously. These models consider bonds of spins 
1 (S - 0, +1) and a crystal lattice field D with phase diagrams exhibiting a tri critical point. 
For instance, the GS model with a magnetic random field (RF) following a gaussian distribution 
has been recently applied to study the effects of the RF on ITs ll22ll . This model is well known 
to present a reentrant first-order phase transition between the spin glass (SG) phase and the 
paramagnetic (PM) one which characterizes an IF. T he p resence of such RF weakens the PM/SG 
first-order transition by lowering the tricritical point [12 3112411 . Most important, it tends to suppress 
the IF. However, such analysis has been done in quite specific conditions: in the limit of strong 
bond disorder and for a continuous distribution of the RF. Thus, a question remains open: are 
those effects found previously still present for other regimes than the strong bond disorder and 
for a discrete distribution of the RF? For both cases one can expect that the phase diagram 
can become richer with the presence of others thermodynamic phases as well as more complex 
boundary phases 1125 1 . 

For this purpose, we studied here the GS model with the presence of a longitudinal magnetic 
RF hi following a bimodal distribution within weak disordered regimes for both the RF and the 
random bonds. Therefore, the model considers a random bond interaction Jij and the RF /i,. The 
interaction Jjj follows a Gaussian distribution with mean Jo/N and variance J^/N, which allows 
one to introduce the bond disorder strength parameter J/Jq in the problem 1 1 8.1 . The disorders are 
then treated by using the replica method, in which the free energy is obtained within the one-step 
replica symmetry breaking (IS-RSB) approximation 112 611 . In the present case, by adjusting the 
disorder strength J/Jq parameter, the strong {J/ Jo » 1) and non-disordered (J/ Jo = 0) regimes 
can be recovered. Most important, the study of J /Jo values corresponding to the weak disordered 
interactions could bring new information concerning mechanisms responsible for the ITs. 

It should also be remarked that the GS model presents the order parameter r (see Eqs. |7]and 
[Tol l related to the diagonal replica elements (r - {S^^S"}). This order parameter can be used to 
obtain the average number of non-magnetic states per site no = ^ - r, which also reflects the 
capability of the sites to interact or not [27.1 . For instance, the limits no = and no = 1 can be 
related to the fully interacting (r = 1 high number of sites with spins S = +1) and the fully no 
interacting (r = high number of sites with 5=0) regimes, respectively. As a consequence, 
there are two types of PM phases with distinct no. This is an important condition for the existence 
of flie IF in the GS model [fll . 

Moreover, the BC model with RF following a bimodal distribution is particularly interesting 
1 2311 . This model does not spontaneously display any IT at mean field level. However, it can 
also present the two types of PM phases with distinct no. It means that one can have a PM phase 
with low entropic content. This fact can be explored to produce an IM in the BC model when an 
artificial entropic advantage of the interacting states is introduced [10]. The IM is characterized 
by a transition from the PM phase at low temperature to the ordered ferromagnetic (F) phase at 
higher temperatures. Nevertheless, the bimodal RF is able to suppress the low entropic PM phase 
1 2311, which could break an essential condition to the occurrence of IT [ 12] in this specific model. 

On the other hand, the regime of weak RB in the GS model also introduces a F phase. How- 
ever, in this is still not completely obvious if an IT can be observed. Thus, one can ask if the 
weak RB regime could be able to produce the necessary entropic advantage of the ordered phase 
to present spontaneously an IT. Particularly this point has not been yet studied in the context of 
the IT problem. In addition, the effects of the bimodal RF could be more evident in the region 
of weak RB disorder. As a consequence, the superposition of the RB (already studied only for 
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strong regimes) and discrete distributed RF (well studied in the BC model) effects on the ITs 
could be more carefully analyzed. To accomplish this task, the effects of both RB and RF in the 
GS model can be traced from the behavior of no ll2lll to the entropic content of thermodynamic 
phases present in the problem. A detailed study of effects caused by different disorders on no can 
be important to clarify the mechanism behind the existence of the IT in the present approach. 

The paper is structured as follows: in Sec. II, the free energy within the IS-RSB scheme 
is found. In Sec. Ill, a detailed discussion of phase diagrams is presented. The last section is 
reserved for the conclusions. 



2. General Formulation 

The model is given by the Hamiltonian 



N N 



(1) 



/=i 



where the spin variables assume the values S - +1,0, the summation (/, j) is over any pair of 
spins and D is the crystal field. The spin-spin coupling Jij and the magnetic fields hi are random 
variables following independent distributions 
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[Jij - Jo/Nf 



nhi) = pSihi -ho) + (l- p)5{hi + ha). 



(2) 
(3) 



Herein the procedure introduced in Refs. 11281 12911 is closely followed. Therefore, the average 
free energy per spin 



/ = - 1 lifiN) j dhidJijP(hi)P(Jij) In Tr exp [-/3HiJij, /i,)] 
can be obtained using the replica method as: 



(4) 



n^O n 



2 ypj] V-' 



2 

'iaa 



Y,^l-P In Tre^' - {I - p) InTre^ | . 



(5) 



where y6 - l/T (T is the temperature), a and y are the replica labels, and are the effective 
Hamiltonians given by 



L — ipj) ^ qayS gS y + ^ ^ QogS 



(6) 
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In Eq. |5]is used the steepest descent method (for ^ oo) to evaluate the fields {m} and [q] that 
are associated with the following correlations 

nia - p < S a >+ +(i - P) < S a >- 

qaa^P< Sl >+ +(1 -p)<Sl (7) 
qafi^ P < SaSp >+ +(1 - p) < SaS/3 >- . 

The symbol < ... >± means the thermodynamic al average over the effective Hamiltonian 
given in Eq. |6] 

The Parisi's IS-RSB is used to parameterize the replica matrix {q} as shown below ll26ll 

q„fi = qiifl(a/x) = l(l3/x) (8) 

qafi^qoifI(a/x)^I(/3/x) (9) 

where /(y) gives the smallest integer which is greater than or equal to y. The replica diagonal 
elements and the magnetization are respectively given by 

qaa — r and - m. (10) 

The parameterization introduced in Eqs. ISTITOlis used in order the sum over the replica labels 
in Eqs. |5]and|6] This procedure results in quadratic terms that are linearised introducing new 
auxiliary fields z and v. Finally, the free energy is given by 

r,f ^-^^^ w i\ 2 2 , 2n , /^Jo 2 



- Jdz In J Dv [K-(z, v)Y , 



(11) 



where 

K^(z, v) = 1 + 2e''cosh/i'=(z, v) (12) 



with JDw^ j ^e-"'/^ (w = z, v), 

y^'-^(r-qi)-l3D, (13) 

and h^(z, v) is 



h^iz, v) = I3J[ y/^z + V?i - gov] +/3{Jom + ho). (14) 

The equations for the order parameters q(),q\,r,m and the block size parameter x are obtained 
from Eq. [TT]by using the saddle point conditions. Other thermodynamic quantities can also be 
derived from Eq. [TT] as for instance, the entropy s - -df/dT. Particularly the present work 
analysis the symmetric case for the RF distribution in which p - 0.5. 




Figure 1 : Phase Diagrams DJJq versus T IJq for ho /Jo = 0.0 and several values of J/ Jo : 0.00, 0.5, 0.75 and 1 .05. The 
full and the dotted lines represent the second and the first-order transitions, respectively. The dashed-dotted lines at low 
T are extrapolations to the T = 0. 



3. Results 

This section displays phase diagrams T/Jo versus D/Jo that are build from numerical solu- 
tions of the order parameters for several set of disorders: J/Jq and h^/Jo. In these, the SG phase 
is characterized by RSB solution (6 = qi - qo > 0) with zero magnetization (m = 0). The ferro- 
magnetic (F) phase occurs when the RS solution {q - q\ - qo) is found with q > and m > 0. 
The mixed ferromagnetic (F') phase is given when m + and 5 > 0. The paramagnetic (PM) 
phase is characterized by m - with RS solution. Particularly the RF induces the correlation 
given by q in the PM phase, but without breaks the replica symmetry ((5 = 0, ^ > with m - 0). 
The behavior of the entropy s and «o fxf\ are also discussed in this section. 

The phase diagrams exhibited in Fig. [T] consider /zq// - 0.0 with several values of the 
strength of bond disorder y/7o : 0.0,0.5,0.75 and 1.05. For example. Fig. [Tta) shows the phase 
diagram without any type of disorder (hoi Jo - 0-0 and JjJo - 0.0). In this case, the continuous 
PM-F phase transition is decreased until a tricritical point as D/Jo increases. After the tricritical 
point, the phase transition becomes first-order where there are more than one PM and F solutions. 
The first-order boundary is then located by comparing the lower free energy of the PM and F (or 
SG) solutions by using the same criteria proposed in Ref. [30J. One important point here is 
that the IS-RSB Ansatz is still unstable against fluctuations in the replica space and the right 
thermodynamics of the SG solutions can be obtained only with full replica symmetry break 
method. However, the IS-RSB procedure can give rehable results around the phase transition 

5 




Figure 2: Behavior of (a) entropy s and (b) no as a function of D/Jq for several values of the strength of bond disorder 
J/ Jo without the presence of RFs. 



(including the reentrance) lHolfTlll . Particularly, the free energy with J/Jq - and ho/Jo - (Eq. 
\TT[ presents the same equation as that one obtained in Ref. II (XI for the Blume- Capel model. In 
fact, the respective phase diagrams found in Refs. jiol 19 1 are recovered in Fig. [Ha) where there 
is no reentrant behavior However, a very weak strength of bond disorder (J/Jq) affects strongly 
the first-order PM-F transition in such way that a reentrant behavior appears. For instance. Fig. 
[Hb) displays a reentrance in the first-order PM-F transition that seems to be intensified as J/ Jo 
increases in Fig. [TJc). Such transition represents an IM, in which the usual entropic relation 
between the ordered F phase and the PM phase is changed (see Fig. |2a)). In addition. Figs. [Tfb) 
andlUc) also indicate the presence of a mixed frustrated F' phase for low T/J, which is obtained 
within the IS-RSB procedure. On the other hand, for the strong disordered regime, the F phase 
is replaced by the SG one (see Fig. [Ttd) for J/Jq - 1.05). In this case there is a continuous PM- 
SG phase transition that decreases until a tricritical point T,c with increasing of D/Jq. Below T,c, 
a reentrant first-order PM-SG phase transition is obtained which indicates an IF. Therefore the 
intensity J/Jq plays an important role in the first-order reentrant transitions. Most importantly, 
the results presented in Fig. [T]suggest that the present model is able to show both the IF and the 
IM. The precise type of IT depends on the strength of the bond disorder 

The mechanism underlying these reentrances can be explained analyzing the entropy s and 
the average number of the non-magnetic sites no- For example. Fig. exhibits s vs D/Jq 
for ho/ Jo - and several values of J/ Jo at low temperature {T/Jo - 0.1). For J/Jq - (case 
without reentrance), the entropy of the ordered F phase is smaller than the PM one at the first- 
order transition. However, the entropy of the ordered phase increases with J/ Jo- In fact, the 
entropy of the PM phase is practically unaffected by increasing J /Jo, while the entropy of the 
ordered phase suffers a large increment (see Fig. |2a) for J/ Jo > 0). Fig. |2](b) shows the behavior 
of riQ for the same values of strength of disorder and temperature as Fig. |2a). This result could 
elucidate why only the entropy of the F phase is modified when the J/ Jo is adjusted. Basically, 
the results indicate that «o of the F phase is increased by J/ Jo- However, no remains unaltered in 
the PM phase that also presents the higher values of no (high number of sites with spin S - 0). 
In other words the low temperature PM phase presents few spin configurations available when 
no — > 1. As a consequence the entropy of this PM phase is very low even for higher J/Jq values. 
It is in contrast with the high temperature PM phase that shows a strong thermal spin fluctuation. 
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Figure 3: Phase diagrams D//o versus T IJq for fto/Zo = 0.4 and Jjjf, : 0.0, 0.5. The full and the dotted lines represent 
the second and the first-order transitions, respectively. The dashed- dotted lines at low T are extrapolations to the T = 0. 
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Figure 4: Behavior of (a) entropy s and (b) no as a function of D/7o for T IJq = 0. 1 , two values of the strength of bond 
disorder J I Jo and two values of hg / Jo : 0.0, 0.4. 



On the other hand, the bimodal RF causes different effects than those ones caused by J/Jo- 
For instance, Fig. |3ja) shows the phase boundary for ho/Jo = 0.4 and J/Jq = 0, where the 
only effect of disorder comes from the bimodal RF. It is important to remark that for J/Jq = 0, a 
coexistence of ordered FM solutions can appear depending on the value of the RF 12311 . However, 
for the purposes of the present work, only the PM/FM border and the effects of the ho/ Jo on it 
are analyzed. In this case, the critical temperature and the tricritical point are displaced for 
lower values of T/Jo and higher values of D/Jo when ho increases (compare Fig. [Ha) with Fig. 
[3a)). Most important, the presence of ho is unable to produce any reentrant behavior Fig. Qb) 
allows one to analyze the effects caused by the superposition of both disorders: ho/ Jo - 0.4 and 
J /Jo = 0.5. It exhibits a reentrance (IM) attributed to the increase of the random bond disorder 
that also introduces a mixed phase at very low temperatures. 

The analyses of the entropy and «o in the presence of RF with J/ Jo = are accomplished in 
Fig|4] Comparing the cases ho/ Jo - 0.0 and ho/ Jo = 0.4 at low temperature {T/Jo - 0.1) in Fig. 
lUa), it is possible to observe that not only the entropy of the F phase increases with ho/ Jo, but 

7 



also the entropy of the PM one. It is important to observe that the entropy of PM phase depends 
on the strength of the RF disorder This increase in s in the PM phase can be one of the reasons 
for the incapabihty of the RF disorders generate IT. The behavior of no in Fig. |3b) can help 
us to understand better this point. For instance, the no is changed by the strength of RF in both 
phases the ordered and the disordered. In particular, the PM phase at lower temperature presents 
a diminishment in no as ho/Jo increases. This is an important difference in relation to the random 
bond disorder As a consequence the RF causes an increase in the entropy of the PM phase. 

4. Conclusion 

In the present work, the Ghatak-Sherrington SG model has been studied with the addition 
of longitudinal magnetic random field (RF) /;,. This study was carried out within the IS-RSB 
scheme using a mean field approximation. This approach can be helpful to elucidate the role of 
different kind of disorders (RB and RF) in the so called Inverse Transitions, melting or freezing. 

According to our results, the crystal field D favors energetically the spin states 5=0, that 
favoring introduces a PM phase for larger D. However, for no disorder (J/ Jq - and ho /Jq = 0), 
the PM entropic content is always larger than that one of the F phase. However, the situation 
changes at very lower temperature when the random bond disorder is considered (J/ Jo > 0). 
For a small strength of bond disorder (0 < J/Jq < 1) the F phase is still found, but now with 
larger entropy as compared with the JjJo - case, while the PM entropy is almost unchanged. 
In other words, it arises an IM for a certain range of D. As the strength of the random bond 
increases {J I Jo > 1) the SG phase replaces the F order and the IF occurs. On the other hand, the 
presence of a RF increases the entropic content not only of the ordered F or SG phases, but also 
of the disordered PM phase. It is important to remark this particular point that the RF causes 
a strong increase of the PM phase entropy. Therefore, it breaks one of the essential conditions 
to the occurrence of IT that is the existence of a PM phase with lower entropy than the ordered 
phase. In fact, the RF tends to couple with the interacting spin states {S - +1) that can be favored 
energetically even for higher values of D. As a consequence, the PM phase at lower temperature 
with smaller entropy is hard to be found. 

To summarize the main difference between the two kinds of disorder studied here in the IT 
context is the effect caused in the PM phase. The RFs compete with D by increasing the entropy 
of the PM phase of lower temperatures, while the random bonds affect only the entropic content 
of the ordered phase favoring the IT. 
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